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Abstract 

In this paper, we propose integrable discretizations of a two-dimensional Hamil- 
tonian system with quartic potentials. Using either the method of separation of 
variables or the method based on bilinear forms, we construct the corresponding 
integrable mappings for the first three among four integrable cases. 
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1 Introduction 



In the last two decades, much attention has been paid to the study of discrete 
integrable systems, which resulted in the discovery of second-order integrable 
mappings, i.e. the Quispel-Robert-Thomson (QRT) mapping [Tf2] and discrete 
Painleve equations [3]. A large number of papers have been devoted to these 
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integrable mappings [Hl5] . In contrast, the examples of higher-dimensional 
integrable mappings are a few pjfTfSPinUfTT] . It is paramount to find and 
study higher- dimensional integrable mappings in order to develop the theory 
of discrete integrable systems. 

One way to achieve above goal is to construct integrable discretizations from 
some known integrable differential equations (DEs). Actually, several pow- 
erful recipes in this aspect have been proposed, such as a method based on 
bilinear form [12] , and a method based on Lax pair P^|14[15|ITB] . Although 
these discretization methods have produced many discrete integrable systems, 
integrable discrete analogues of some many well-known integrable higher- 
dimensional (coupled) Hamiltonian systems are still missing. 

In this Letter, we construct integrable discretizations of a two-dimensional 
Hamiltonian system with a quartic potential, whose Hamiltonian can be ex- 
pressed as 



H = ~ (pj + p|) + Aql + Bq\ + aq\ + (5q\ + bq\q\ . (1) 
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Thus, the equations of motion are given by 

<?2 = g = P2 , 

(2) 

Pi = "If = - 2 ^i - 2fc(2ag? + 5q 2 ) , 
P2 = -ff = ~2Bq 2 - 2q 2 {2(3q 2 + Sqf) , 

or 

q\ = -2A qi - 2q x (2aq 2 + 5q 2 2 ) , 

(3) 

g 2 = -2Bq 2 - 2q 2 {2f3q 2 + 6qf) . 

The above Hamiltonian system is known to be integrable for the following four 
cases [T7|18|19j : 

(I) P = a, 5 = 2a, a, A, B arbitrary, 

(II) /3 — a, 8 — 6a, A = B, a, A arbitrary , 

(III) (3 = 16a, 5 = 12a, B = 44, a, A arbitrary, 

(IV) (3 = 8a, 5 = 6a, B = 44, a, A arbitrary , 

with the corresponding second integrals listed as follows: 

(I) I( qi ,q 2 , Pl ,p 2 ) = (qiP2-q2Pi) 2 + ^{p 2 1 + 2qj[A + a(q 2 1 +q 2 2 )]}, 

(II) I{qi,q 2 ,Pi,p 2 ) = P!p 2 + 2q 1 q 2 [A + 2a{q\ + ql)} , 

(III) I{qi,q 2 ,Pi,p 2 ) = (qm- Q2Pi)pi + 2q 2 q 2 [A + a{q 2 + 2q 2 )] , 

(IV) I(q 1 ,q 2 ,p 1 ,p 2 ) = Aaq 2 { qi p 2 - 2q 2Vl f + {p\ + 2q 2 [A + a(q 2 + 2q 2 )]} 2 . 
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The cases (I), (II) and (III) are separable in ellipsoidal, Cartesian and paraboidal 
coordinates, respectively, while the case (IV) being separable in the general 
sense [20|21j22j. Further, the case (I) is equivalent to a traveling wave reduc- 
tion of the integrable coupled nonlinear Schrodinger (cNLS) equation (Man- 
akov system), and it also corresponds to the Gamier system of two particles 
[23p4]- 



The present paper is organized as follows: In section 2, we present two dis- 
cretizations of the case (I) based on bilinear form. Then, we show the integra- 
bility of these discretizations by deriving discrete versions of two first integrals 
for the original Hamiltonian system (CTJ. It should be noted that what we ob- 
tain here are alternative forms of some already known results [25 2 6]fl6] . In 
section 3, we derive an integrable discretization of the case (II) and show its 
integrals. The key of the discretization is the separation of variables. In section 
4, we construct an integrable discretization of the case (III) by virtue of the 
separation of variable. The approach used in this case is similar to the one 
used for the discretization of the Kepler problem [27J. We give some remarks 
and further topics in section 5. 
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2 The discretizations of the case (I) and Integrals 



2.1 The discretizations of the case (I) 



In this section, we consider the integrable discretization of the case (I), i.e. 
(3 = a, 5 = let (a, A and B: arbitrary). Upon the dependent variable trans- 
formation 



ft = j, i = l,2, 



eq. (EJ) can be casted into the bilinear forms 



r% gi .f + 2Afgi = 0, 
-D 2 t f-f + Aagl + Aagl = U, 
D 2 t g 2 -f + 2Bfg 2 = 0, 
-D 2 t f-f + Aagl + Aagi = 0. 



(4) 
(5) 
(6) 
(7) 



Based upon the discetization of above bilinear forms, we can obtain two inte- 
grable discretizations as follows. 

Discretization (a) 

An integrable discretization of eq.® for the case (I) is 



Ql,n+1 ~ 2qi,n + Ql,n-1 



h 2 



-2Aq 1>r 



-2aq hn 

Q2,n+1 — 2q 2: n + <72,n-l 



I -Ah 2 

9l,n(?l,ti+l + Ql,n-l) + i ^?2,n(?2,n+l + ?2,n-l; 



1-Bh 2 



h 2 



-2Bq 



2,n 



-2aq 2 , n 



Q2,n{Q2,n+l + Q2,n- 



1-Bh 2 
I- Ah 2 



Ql,n (ll,n+l + ?l,n-l) 



(8) 



(9) 
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Upon the dependent variable transformation 



ft,n = % i = 1,2, (10) 



eqs.dS])-© are transformed to 



01,n+l/n-l + 9l,n-lfn+l = ^l9l ,nfn > (H) 

1 _ Ah 2 

(2 - 2Ah 2 )f n+1 f n _ 1 - 2a/i 2 r 1 ^ jB - 2c* T — ^pT^ = TJ 2 , (12) 

92,n+lfn-l + 92,n-lfn+l = ^292,nfn , (13) 

1 - £?/i 2 

(2 - 2Bh 2 )f n+1 f n ^ - 2ah 2 Y 2 g\ n - 2a — ^ h 2 T l9 l n = T 2 f 2 , (14) 



which are discretization of bilinear forms (j4])-(J7j). 

Multiplying eq.® by (gi,„+i - gi,„-i)/2, eq.© by (g 2 ,n+i - &,n-i)/2, and 
summing them up, we obtain an integral of motion 



#™ = ^Pi,n + \v\n + Aq 1<n+1 q hn + Bq 2 , n +iq2,n + aqln+illn 
2 2 (I- Ah 2 1-Bh 2 \ 

+ a( l2,n+l ( l2,n + " I j _ ^ 2 + j _ ^ 2 I 9l,n+l?l,n?2,n+l92,n , (15) 



which is a discrete analogue of the Hamiltonian ([!]) where 

<?l,n+l — <Zl,n _ Q2,n+1 ~ Q2,n 

Pl,n- ~ h , P2,n~ ~ h ■ 

Discretization (b) 

Another discretization of eq.([3J), which was firstly proposed by Suris [25lfT6] 
is of the form 



q ln+1 - 2q ln + gi 

— = -2Aq hn - 2a(q 1>n + q 2> „){qi, n +i + ?i,n-i) , (16) 
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H2,n+l — ^2,n T q2,n-l OD r> / 2 , 2 \ / , \ /-.^ 

— = -2Bq %n - 2a{q l n + g 2jn )(g 2 ,fH-i + Q2, n -i) ■ ( 17 ) 

Upon the same transformation (ITU1) . eqs. (ITBl) - (fT7j) can be casted into discrete 
bilinear forms 

(18) 

fl , (19) 

(20) 

/n- (21) 



91, n+lfn-l + 9l,n-lfn+l — ^l9l,nfn i 

1 fn+ifn-i ~ 2ah 2 g\ - 2ah 2 g\ 
1 1 

92, n+lfn-l + 92,n-lfn+l = ^292,nfn , 

p" 5 ^ fn+ifn-i ~ 2ah 2 g\ n - 2ah 2 g\ n 



Similarly, it can be shown that discretization (b) admits an integral of motion 

i / 2 2,2 2,2 2,22 \ /'OOA 

+ a Wl,ri+l9l I n + l2,n+ll2,n ' Q2,n+lQl,n ' Q2,n1l,n+l) J 

which is a discrete analogue of the Hamiltonian (pQ). 



T/ie integrals of discretizations of case (I) 



It was shown by Suris [251TT6] that the discretization (b) admits two integrals 



h,n = qln+i + 5?,„ ~ 2 (1 ~ Ah 2 )q hn+1 q hn + 2ah 2 q 2 ln+l q 2 ln 
2a(l-Ah 2 ) 2 2 2 , 2 2 , 

(1 - Ah 2 ) 2 - (1 - Bh 2 ) 2 ^ 1 > n + l92 '« %n%n+l) 
Aa(l - Ah 2 ){l - Bh 2 ) 
_ (1 - Ah 2 ) 2 - (1 - Bh 2 ) 2 ^ l ' n+iq2 ' n+xqx,nq2 ' n ' 
hn = (&,n+i + (&,n ~ 2 (1 ~ Bh 2 )q 2)n+1 q 2 ^ n + 2ah 2 q\ n+l q\ n 
2a(l-Bh 2 ) 2 2 2,2 2 \ 

(1 — Bh 2 ) 2 - (1 - Ah 2 ) 2 ' n+iqi ' n q2 >n q hn+i) 



(23) 
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Note that the sum of above two first integrals corresponds to the discrete 
analogue of Hamiltonian ( 1221) , the difference corresponds to the second integral 
of motion. 

Next, we show that discretizations (a) and (b) share the same integrals. Setting 
I\ = 1 — Ah 2 and = 1 — Bh 2 in the bilinear equations (|TT|) - (fl3|) . we have 

(1 - Ah 2 )q lin (q 2t n+i + <?2,n-i) = (1 - Bh 2 )q 2 , n {qi,n+i + <7i,n-i) • (25) 

Employing ( f25|) . it is easily shown that eqs.®-© are equivalent to eqs.lfTB"])- 
(fl71) . Thus, discretization (a) has the same integrals as discretization (b). 

Special case (A = B): 

Under this case, the discretization is 

gi jn+ i + gi jn _! - 2(1 - Ah 2 )q ltn + 2ah 2 (ql n + q% in )(qi, n +i + qx,n-i) = 

(26) 

<?2,n+i + q%n-\ ~ 2(1 - Ah 2 )q 2) n + 2ah 2 {ql n + ql n )(q 2) n+i + q2,n-i) = . 

(27) 

Multiplying eq.^M) by {q 2 , n +i + <?2,n-i), eq.([2ZD by {qi, n+ i + qi, n -i), and sub- 
tracting them, we arrive at 

-2(1 - Ah 2 )q hn (q 2tn+1 + g 2 , n -i) + 2(1 - Ah 2 )q 2tU (q hn+1 + gi, n _i) = , 
from which, we have the following second integral 

I = 9l,n92,n+l — 9l,n+l92,n • 



It is known that the case (I) can be obtained by a stationary reduction of a 
Manakov system. This fact remains true for the discretization, i.e. the inte- 
grable discretization of the case (I) can be obtained from a semi-discrete cNLS 
equation: 



-2Aq ljTl - (2a\q hn \ 2 + 2a\q 2tn \ 2 )(q 1 

,n+l + Ql ,n- 



• d _ I 92,n+l— 2g2,n+92,n-l 

1^92,n H p 



f28) 



-2Bq 2j7l - (2a\q hn \ 2 + 2a\q 2>n \ 2 ){q 2 ,n+l + ?2 



n-1; 



by a stationary reduction. The semi-discrete cNLS equation ( 1281) was studied 
by Ohta and iV-soliton solutions were given in the form of Pfaffian [26J. 

3 Integrable discretizations of case (II) 

Let us consider how to build up integrable discretization of the case (II) {(3 = 
a, 5 = 6a, A = B), where the equations of motion are given by 

q\ = —2Aqi - 2qi(2aq\ + Qaqj) , 

q 2 = -2Aq 2 - 2q 2 {2aq\ + Qaqj) . (29) 

First, we comment that the method of discretization based on bilinear form is 
not able to construct integrable discretization of the case (II). This is a result 
of breaking solution structure of the case (II) by the transformation to the 
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bilinear forms. 

Fortunately, an integrable discretization can be constructed by noting that 
eqs. (1290 are separable in Cartesian coordinates. It is known that eqs. f!29|) can 
be decoupled to two anharmonic oscillators by a simple transformation: u\ = 
q 1 + q 2 , u 2 = qi - q 2 : 
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iii — —2Au% — Aau 1 



u 2 = —2Au 2 — Aau\ . 



An integrable discretization 



Wl n +l — 2U\ n + Mi „_i a 2 U l,n+1 + u l,n-l 

— = -2Au 1}U - Aau l n , 

U 2 n+1 — 2u 2n + U 2n -\ 2 u 2,n+l + U 2 n^\ 

— = -2/4w 2 ,n - ™ 2 ,n ^ ' ( ) 



is known, then we can easily give two integrals which are discrete Hamiltonians 
of ([30]) as follows: 



Hi = -p^, n + Au hn u 1<n+1 + au\ n u\ tn+l , 

#2 = ^4,„ + Au 2)n u 2in+1 + au\ n u\ n+1 , (31) 



Substituting U\ >n = qi tn +q 2 , n and u 2 ,n = qi, n —q2,n into two discrete anharmonic 
oscillators (13"U|) . we have the following discretization of the case (II) (symmetric 
discretization): 

<?l,n+l — 2gi >n + qi yn -i 



h 2 

-2Aqx tn - (aq hn + aq 2 ) 2aq ljU q 2 



2 2 \ + Ql,n-1 <?2,n+l + <?2,n-l 
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Q2,n+1 — 2^2,™ + l2,n-\ 



h 2 

-2Aq 2) n - («gi n + aq 2 ) 2ag 2 ,n9i 



J 2 \ l2,n+l + 92,71-1 9l,n+l + 9l,n-l 



(32) 



Similarly, we can obtain an asymmetric discretization of the case (II) by sub- 
stituting ui >n = 9i, n + 92,n an d M2,n-i = Qi,n ~ Q2,n into two discrete anharmonic 
oscillators fl30l) 

9i,n+i — 2gi >ra + qi t n-i 

nA i 2 2 \ 9l,7i+l + 9l,n-l ?2,n + ?2,n-2 

= -2Agi, n - (ag 1>n + ag 2 ,n-i) ^ 2agi jn g 2 ,n-i ^ , 

Q2,n+1 — 2^2,™ + l2,n-\ 
h 2 

nA / 2 2 \ 92,n+l + ?2,n-l 91,71+2 + 91,71 

= -2Ag 2 ,n - (ag 1>n+1 + ag 2 ,J ^ 2ag 2 ,n9i,n+i g ' 

(33) 



It can be easily shown that (ifi + H 2 )/2 and (i/i — H 2 )/2 in (13T|) are discrete 
versions of two integrals for the original Hamiltonian system (pD) 



J l = \v\n + T^U + ^(9l,n9l,n+l + 92,n92,™+l) 

+«(?! 9l,7l+l + 9 2 ,n9 2 ,71+1 ,71+1 ) J 

^2 = Pl,nP2,n + A(gi jn g 2 ,n+1 + 92,n9l,n+l) 

92,n+l + 92,n9l,7i+l )(<& 9l,7J+l + Q2,nl2,n+1 ) , (34) 



by the dependent variable transformation n l n = g l ra + g 2 ,n and u 2jn _i = 

91,71 — 92,71- 

In Ref . [B] , the more general form of integrable discretization of the anharmonic 
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oscillator 

d 2 x 



, + ax + 6x 3 = 
at 2 



was presented. The form of the discretization of the anharmonic oscillator is 



u n+l 2li n + Un-1 i o /I r / i M 
— 2 h 2A[c U U n + Ci 2 {U n+ i + U n -1 j J 

+4a[c 2 iM n+ iM„,M„_i + C 2 2U 2 n {u n +l + «n-l)] = ■ 



Using this discretization and the transformation -ui jn = gi jn + g 2) „, W2,n = 
<?i,n — <?2,n, w e can make a more general integrable symmetric discretization of 
the case (II), 



ji = -2-4 (cngi jn + c 12 {qi, n +i + <7i,n-i)) 

-4ac 22 (gi in + 2aql n ){q ljn +x + gi,„-i) - 8ac 22 g lin (g 2in (g 2>n+ i + ?2,n-i)) 

— 2ac 21 [g2,n92,n+l W.n+l + <Zl,7i-l) + 9l,n9l,n+l9l,n-l + <?2,n<?2,n+l<?2,n-l] , 
#2,71+1 — 2g 2 ,„ + 92,71-1 n A / ( , \ \ 

75 = -2 A (ciig 2 ,n + ci 2 (g 2 , n+ i + g 2 , n -i)) 

-4ac 22 (g 2n + 2aql n ){q 2 , n +i + <?2,n-i) - 8ac 22 g 2jn (gi in (gi in+ i + gi, n -i)) 

— 2«C 2 i [gi, n 9l,n+l (92,n+l + ?2,7i-l) + 9l,n9l,n+l9l,n-l + 92,n92,n+l92,n-l] • 

(35) 



The constants en, C\ 2 , c 2 \ and c 22 are given in Ref. [6j, i.e. 

cn=7ii, ci 2 = 7i2, C21 = (1 - 07u/2)7 2 i , c 22 = (1 + a7i 2 )7 22 /2 

while 7u,7i2,72i,722 satisfying the relations 7n + 27 i2 = 1, 721 + 722 = 1- 
From the result in Ref.|53, we can construct integrals of eq. (13"5"|) : 

h = H 1 + H 2 , 
I 2 = Hi — H 2 , 
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where 



Hi = (Cu[ci 2 ((gi,n + <?2,n) 2 + (?l,n-l + <?2,n-l) 2 ) - Cn(gi, n + g 2 ,n) + <?2,n-l)] 

+2a(cnc 22 + cnc 2 i)(gi, n + g 2 ,„) 2 (gi,„_i + <?2,n-i) 2 ) 
/(cii[cn + 2ac2i((gi,„ + g 2 , n ) 2 + (?i,n-i + <?2,n-i) 2 )] 
+4a 2 c 21 (gi, n + g 2 ,„) 2 (gi,„_i + g2,n-i) 2 ) , 

#2 = (Cu[ci 2 ((gi,n - ?2,n) 2 + (?l,n-l ~ <?2,n-l) 2 ) ~ Cn(<7i >n - g 2 ,fi) (ft.n-l ~ ?2,n-l)] 
+2a(cnc 22 + cnc 2 i)(gi ) (tfl.n-l - ?2,n-l) ) 

/(cn[cu + 2ac2i((gi,„ - g 2 ,„) 2 + - <?2,n-i) 2 )] 

+4a 2 c 21 (gi in - g 2 ,„) 2 (gi,„_i - <?2,n-i) 2 ) • 



4 Integrable discretization of case (III) 



Let us consider the case (III) (f3 = 16a, 5 = 12a, B = AA), whose Hamiltonian 
is 



H = \ (pi + Pi) + Aql + AAq 2 2 + aq\ + 16aq\ + 12aq\q 2 2 . (36) 



The equations of motion are therefore given by 
q\ = -2A qi - 2q x {2aq\ + 12aq\) , 

(37) 

q 2 = SAq 2 - 2q 2 (32aq\ + 12ag 2 ) . 

Eqs. (I37p are also separable. Upon the canonical transformation called the Levi- 
Civita transformation 



q 1=Ul u 2 , q 2 = -( Ul -u 2 ), pi = — — , p 2 = — — 2 ,(3S 

2 u( + u 2 u{ + u 2 



13 



eq. (l36|) is converted to 



H(p ul ,p u2 ,u 1 ,u 2 ) = 2 ^ 2 (\{p 2 ul + pis) + A(ut + u%) + a(u{° + v™)) .(39) 



By taking the level set A of Hamiltonian, i.e. 



H = A 



and introducing a new variable s defined by 
dt = {u\ + u\)ds 



> 



we obtain 



^=Puk, ^r- = ~2\u k + 6Aul + Waul , k = l,2, (40) 
as as 



(41) 



or 



<ffi = -2Awi + QAul + Waul , 
^ = -2\u 2 + QAu\ + I0au 9 2 ■ 

The discretization of eq. (|40|) is given as follows: 

Uk,n+1 Uk,n Puk,n+1 Puk,n 

s (n+l) ~~ S (n) 2 

3 3 5 5 

Puk,n+1 ~ Puk,n ~ , , n , 3 U k,n+1 ~~ U k,n-1 , 5 U k,n+1 ~ U k,n-1 

— : = -2Xu k , n + 2Aufc n — ■ ■ h 2aw° n — ! ■ , 

s (n+l) ~~ S( n ) ' Ufc )T i+l — Uk, n -1 ' Uk, n +1 ~ Uk,n-1 

k=l,2, (42) 
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where 



Pu k ,n = — " 

s (n+l) — S( n ) 



It is shown that the integrals are 



7- ^ 2 ^ /t 3 3 5 5 

-'fc = -^Puk,n + ~^ u k,n+lUk,n + ^ u k,n+l U k,n + u k,n+l u k,n 



Introducing a discrete variable t n as in 



eq. (|42p can be rewritten as 

Uk,n+1 ^k,n Puk,n+1 Puk,n 

t(n+l) - t( n ) 2(«f n + u\ n ) 

Puk,n+1 Puk,n 



t(n+l) — t(n) 



' ^ \ i o /t 3 M £i,n+1 U k,n-1 . 5 

-2Xu kin + 2Au kn — ! ! h 2au{ 



{u\ n + u\ n ) \ "'"U^n+l - Uk, n -1 

k = 1,2. 



Since 

1 /l 



we have 

Uk,n+1 ^k.n Puk,n+1 Puk,: 



t(n+l)—t(n) 2( u l,n + u 2,n) 
Puk,n+1 Puk,n 

t(n+l) ~ t(n) 
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X 




.71. 



(45) 



Using the Levi-Civita transformation (I38I) . we can compute values of q±, q2,Pi,P2- 
5 Conclusions 

In this paper we have searched for integrable discretization of a two-dimensional 
Hamiltonian system with a quartic potential. Among the four integrable cases, 
we have succeeded in constructing integrable discretizations for the first three 
cases, i.e., (I), (II) and (III). The key of integrable discretization is the sepa- 
ration of variables. 

Since our method of discretization using separation of variables can be ap- 
plied to many higher-dimensional Hamiltonian systems, it is possible to obtain 
various integrable discretizations of higher- dimensional Hamiltonian systems. 
These aspects will be discussed in the forthcoming paper. 
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